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Abstract 



o 

Q 

An adaptive optimization technique to improve precision of quantum ho- 
modyne tomography is presented. The method is based on the existence of 
so-called null functions, which have zero average for arbitrary state of radi- 
' ation. Addition of null functions to the tomographic kernels does not affect 

their mean values, but changes statistical errors, which can then be reduced 
■ by an optimization method that " adapts" kernels to homodyne data. Applica- 

CN , tions to tomography of the density matrix and other relevant field-observables 

\ are studied in detail. 

Q\ 
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I. INTRODUCTION 



The possibility of measuring the quantum state of radiation has been received an increas- 
ing interest in the last years |l|-|3| , as it opens perspectives for a new kind of experiments in 
quantum optics, with the possibility of measuring photon correlations on a sub-picosecond 
time-scale ||], characterizing squeezing properties ||, photon statistics in parametric fluo- 
rescence @, quantum correlations in down- conversion 0, nonclassicality of states ||, and 
measuring Hamiltonians of nonlinear optical devices || . Among the many state reconstruc- 
tion techniques suggested in the literature [|T0|— pT9H , quantum homodyne tomography (QHT) 
p~T| ^T~3| , p~Hj] of radiation field have been received much attention being the only method 
which has been implemented in quantum optical experiments and recently being 



extended to estimation of the expectation value of any operator of the field ]TB[ , which makes 
the method the first universal detectors for radiation. 

On one hand, QHT takes advantage of amplification from the local oscillator in the 
homodyne detector, avoiding the need of single-photon resolving photodetectors, hence with 
the possibility of achieving very high quantum efficiency using photodiodes 0. On the other 
hand, the method of QHT is very efficient and statistically reliable, and can be implemented 
on-line with the experiment. 

In principle, a precise knowledge of the density matrix would require an infinite number of 
measurements on identical preparations of radiation. However, in real experiments one has 
at disposal only a finite number of data, and thus statistical analysis and errors estimation 
are needed. The purpose of this paper is to analyze the possibility of improving the current 
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QHT technique, in order to minimize statistical errors. We will present a new method that 
"adapts" the tomographic estimators to a given finite set of data, improving the precision 
of the tomographic measurement. 

Quantum tomography of a single-mode radiation field consists of a set of repeated mea- 
surements of the field-quadrature x^ = \{ae^ + a^e 1 ^) at different values of the reference 
phase (f). The expectation value of a generic operator can be obtained by averaging a suitable 
kernel function R[0](x,<f)) as follows |]TB 



(6) = Tr \ed}= rdxp v {xA)Rn[0]{xA) , (1) 

L } JO 11 J-oo 

where p v (x, 0) denotes the probability distribution of the outcomes x for the quadrature x$ 
with quantum efficiency rj, and R v [0}(x,4>) is given by 

R v [6](x,<P) = - -! drexp ^-p-r Tr \d cos [y/¥(x - x+)] } . (2) 

4 Jo L °v J ^ J 

In the following we will focus attention only on the case r\ = 1, and we will drop the subscript 
7] in the notation. As it will appear from the following, the method works equally well also 
for nonunit quantum efficiency, and a detailed numerical analysis versus 77 will be given 
elsewhere. On the basis of identity (JT|) , it follows that the ensemble average (O) can be 
experimentally obtained by averaging R[0](x,(j)) over the set of homodyne data, namely 



1 N 

(0) = R[0] = -J2 R {°\^iAi) , (3) 



N being the total number of measurements of the sample. The statistical error of the 
tomographic measurement in Eq. (^j can be easily evaluated provided that the corresponding 
kernel function satisfies the hypothesis of the central limit theorem, which assures that the 
partial average over a block of data is Gaussian distributed around the global average over all 
data. In this case, the error is evaluated by dividing the ensemble of data into subensembles, 
and calculating the r.m.s. deviation of each subensemble mean value with respect to the 
global average. The estimated value of such a confidence interval is given by 

1/2 



50 = - 7 =\AR*[0]\ . (1) 



where Ai? 2 [0] is the variance of the kernel over the tomographic probability 

AR 2 [0} = J°°dx p(x, 0) R 2 [6] (x, 0) - I J°°dx p(x, 0) R{6] (x, 0) J . (5) 

Following this scheme, the tomographic precision in determining matrix elements of the 
density operator g has been discussed in with asymptotic estimations in Ref. 

|2^] , whereas relevant observables O have been analyzed in [23], also in comparison with the 
corresponding ideal measurement. 



2 



The crucial point of the method presented in this paper is that the tomographic kernel 
R[0}(x, 0) is not unique, since a large class of null functions p4,25| F(x, <p) exists that have 
zero tomographic average for arbitrary state, namely 



Jo * J-oo 

Therefore, addition of null functions to a generic kernel gives a new kernel with the same 
tomographic average, hence equivalent for the estimation of the same ensemble average (O). 
On the other hand, adding null functions would modify the kernel variance, whence the 
statistical error over data. The adaptive tomography method thus consists in optimizing 
kernel in the equivalence class, in order to minimize the statistical errors. 

The paper is structured as follows. In Section || we introduce the classes of null functions 
that will be used in the paper, and describe the adaptive optimization method in detail. In 
Section |TTl| we apply the adaptive method to the tomography of the density matrix in the 
photon number representation. In Section |I^ we analyze the improvement of precision in 
tomographic measurement of some relevant field-observables. Section |V| briefly describes the 
effects of systematic errors on the effectiveness of the method. Finally, Section [VT| closes the 
paper by summarizing the main results. 



II. ADAPTIVE TOMOGRAPHY 

The following functions have vanishing tomographic expectation (g) 

G+ {x,<f)) = e i{l+n)2,i> g + {xe i(t> ) G~ (x,<j>) = e - i{l+n ^ (xe"^) . (7) 

In Eq. ([7|) n > and g±{z) are analytic functions of z. The set Q of null functions defined in 
Eqs. (^) forms a vector space over C, and each class Q ± = {G^} separately is closed under 
multiplication (without inverse). 

In order to prove vanishing expectation (|6|) for G^(x, <fi) we consider the Taylor expansion 
of functions g±(xe' 1 ^) 

oo 

g±{xe i +) = 4x k e ±ik<t> , (8) 

k=0 

which allows to write 

G|= rdxp{x,(i ) )e ± ^ +n ^g ± {xe^) = Yct e ± ^ k+2+2n ^ (zj) , (9) 

JO * J-oo Jo * 

where (■ ■ ■) denotes the usual ensemble average. Using the Wilcox decomposition formula 
|26| one can write 

[[fc/2]] k-2p f . fc _2p_ a » 

(Ak\ - V V ^ a I J(2p+2s-k)<f> / 1n x 

* 2fc ^ ho 2p p lsl ( k - 2 P~ S V- ' 

where [[x]} denotes the integer part of x. Eq. ([U]) together with the identity 
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q even 



e *i<P = J 1 q = , (11) 
1 ^ q odd 



prove that 



hence 



r# e±i ( fc+ 2 + 2n)0 ^ = ) n > , A; > , (12) 
/ ^ / dxp(x,<P) G%(x,<j>) =0 , n>0, (13) 

JO 7T ./-oo 

namely Cr„ (x, 0) are null functions for n > 0. 

In the following, we will focus attention on three particular sets of null functions. The type-I 
null functions are obtained from Eq. ([?]) by choosing n = and g{xe l< ^) = x k e lk ^ for a given 
fc, and will be denoted by F[ (x, 0), namely 

F fc 7 (x, 0) = xV (fc+2) * fc = 0, 1, ... . (14) 

The type-II null-functions correspond to the simple choice g(xe l( ^) = 1, i. e. 

F„ /7 (0) = e i(1+n)2< ^ n = 0,l,... . (15) 

Finally, the type-Ill null functions are a kind of intermediate choice between type I and type 
II classes, and are defined as follows 

F{ n {x, 0) = x k ^ e *m+*+Mtl)4> I = 0, 1, ... , (16) 

where k[l] and n[l] are given in Table |. In the following we will use the notation Fk(x,<p), 
dropping the type index I- 1 1 1 , when the identity under consideration holds for all three types. 
Let us consider a generic real kernel R[0](x,<f)). By adding M null functions keeping the 
kernel as real, we have a new kernel X[0](x, 0) 

M-l M-l 

K[6](x,<P) = R[6](x,<f>) + VkF k {x,4>) + Y, »l*kM) , (17) 

k=0 k=0 

where F k (x, 0) G Q + , F£(x, 0) G Q~ , and \i k are complex coefficients. By definition we have 
K[0] = R[0], whereas the variance of the new kernel K[0](x, 0) is given by 

'M-l M-l M-l 



AK"[6} = AR2[6] +2\J2 WiFuF? + £ ^R[0]F k + £ fi* k R[0]F* . (18) 



,fcj=0 fc=0 fc=0 



In deriving the above formula we use the fact that both Q + and Q~ are closed under multi- 
plication. 

The variance of the modified kernel function in Eq. ( |I8|) can be minimized with respect 
to the coefficients [i k , leading to the linear set of equations 
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£ t lnF k F? = -R[d]F;. (19) 
i 

It is convenient to rewrite the optimization equations ( |i~9l) in matrix form as follows 

A/i — b. (20) 
where A is the Hermitian M x M matrix 

A kl = F\Ff= [ ^ f dxp(x,(f>) F k (x,<f>)F^x,<f>) , (21) 

Jo 7T J-oo 

and b is the complex vector 

b k = -R[0]F* = - f ^ / ^^^^(^^^(x^). (22) 

JO ^ J-oo 

Notice that the vector b depends on both the kernel R[0] and the state g under examination, 

whereas the matrix A depends on the state only. 

By substituting Eq. (|H]) in Eq. (|H|) and inverting Eq. ( |2"0"D we obtain 

M-l A/-1 

A 2 [6] = AR?[6] - AK^O] = 2 »k A M ^ = 2"£ h * 6? > , (23) 

k,l=0 k,l=0 

which expresses the variance decrease in terms of A and b. 

Let us summarize the optimization procedure for the kernel R[0](x,<f>). After collecting 
an ensemble of iV tomographic data, the quantities A and b are evaluated as tomographic 
experimental averages. Then, by solving the linear system ( p0| ) one obtains the coefficients 
fik which are used to build the optimized kernel K[0](x, 4>). At this point, the same data set 
is used to average K[0](x, (ft) and, upon dividing the set into subensembles, the experimental 
error is evaluated, whose square now is reduced by the quantity A 2 [0]/N. 
The actual precision improvement of the tomographic measurement depends both on the 
state under examination (which affects both b and A) and on the operator O, whose kernel 
enters only in the expression of b. An explicit expression for Am can be obtained by means 
of Eq.(|H]), and generally depends on the type of null function that are involved. For type-II 
null functions it reduces to the identity matrix, independently on the state 

A-m = 5m type-II null functions , (24) 

8m denoting Kronecker delta. For type-I null functions one has 

< = W" g 2, P Hi- P y.(k- P y. type - 1 nul1 funct,ons ' (25) 

The explicit expression for coherent and Fock states is 

A H = a*- ^ fc *^ ! 2- k - 2l L k l -\-2\a\ 2 ) coherent state \a) (k > I) (26) 

<r\k— n+1 poo 

A!* = 5 M / dy e"^ y 2k H 2 (y) Fock state \n) , (27) 

ntyir Jo 

where H n (x) denotes Hermite polynomials. Notice that for Fock states the matrix is diagonal 
(which is true also for type-II and type-Ill null functions). 
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III. ADAPTIVE TOMOGRAPHY OF THE DENSITY MATRIX 



In this section we apply the adaptive method to the tomographic measurement of the 
density matrix in the photon number representation. We evaluate the variance reduction 
A 2 [|n)(m|] in Eq. (^) for O = \n)(m\ corresponding to the tomographic measurement of 
the matrix elements g nm = (m\g\n). We consider the different types of null functions, and 
calculate A 2 [|n)(m|] versus the number M of added null functions, for either coherent states, 
squeezed vacuum, Fock states, and the " Schrodinger-cat" like superposition of coherent 
states given by 



1 



2W1 + exp(-2|a| 



| a) + | — a) 



(28) 



In order to see the new adaptive method at work Monte Carlo simulated experiments are 
presented. 

Tomographic kernels for the matrix elements in the Fock basis have been firstly presented in 
Ref. |fn^| , with extension to non unit quantum efficiency in Ref. [13[|, and factorization iden- 



tities for the kernel in Ref. |p7|| . However, none of the above methods allows for an explicit 
analytical evaluation of the vector b in Eq. (|22|). For this reason, we compute A 2 [|n)(m|] 



numerically presenting results in terms of the relative variance reduction 7, defined as follows 



AK 2 [0] A 2 [0] 



7 = 1 

AR 2 [0} AR 2 [0] 

A complete removal of fluctuations would correspond to 7 = 1. 



(29) 



A. Coherent States 

The adaptive method leads to a significant error reduction for detection of matrix ele- 
ments (m\g\n) of coherent states. Our results indicates that type-I null functions are the 
most effective, and that the larger is the amplitude a of the coherent state, the larger the 
noise reduction. In Fig. [l| numerical results are presented for diagonal elements (n\g\n) for 
intensity |a| 2 = 5. In Fig. |I|(a) the noise reduction 7 is given versus the number M of added 
type-I null functions. One can see that the noise reduction 7 saturates for large M, and 
better levels 7 of reduction are achieved for smaller n. In Fig. |I](b) the noise reduction is 
reported versus n for M = 30. In Fig. |2| we report the results from a Monte Carlo experi- 
ment for I o; 1 2 = 3, with optimization performed with M = 6 null functions. The reduction 
of statistical errors for low values of n is evident. 

The noise reduction for the off-diagonal matrix elements behaves similarly to the diagonal 
ones, being more effective for low indices. In Fig. [3] the noise reduction 7 versus n and m of 
the matrix element (m\g\n) is plotted for a coherent state with \a\ 2 = 10, and for the three 
types of null functions. The type-I null functions are generally more effective, though not 
uniformly over all indices n and m. 
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B. Squeezed states and Schrodinger cat states 



Results for squeezed states and "cat" superposition of coherent states are presented in 
the same subsection, since they behave similarly. This is due to the fact that both states 
have phase- dependent features, that reflect in a similar odd-even oscillation in the photon 
number probability distribution. In Figs. |] and [5] the noise reduction for both cases is 
plotted for the three types of null functions, for M = 10. From the plots it is apparent 
that type-II null functions are now the most effective ones, especially for off-diagonal matrix 
elements, though the same level of noise reduction for low n and m can also be obtained 
using type-I and type-Ill null functions. In Fig. |^ results from a Monte Carlo simulated 
adaptive tomography on a squeezed vacuum are reported for (n) = 4 and M = 10. Matrix 
elements before and after optimization can be compared, showing the error reduction at 
work. 



C. Fock states 



For Fock states the matrix A is diagonal for all types of null functions, and therefore the 
optimization procedure just consists of the evaluation of the vector b. The kernels for the 
matrix elements have the form K[\n) (m\)(x, <fi) = f n ,m{ x ) exp(z(n — m)<f)), where f n , m ( x ) has 
the parity of n — m P,27|. This fact, together with the integral ([11]) makes straightforward 



to show that 



bl = hi 1 = b{ n = Vk , (30) 



namely no improvement should be expected for the precision of quantum tomography on 
Fock states. 



IV. ADAPTIVE TOMOGRAPHIC MEASUREMENTS OF OBSERVABLES 

The tomographic estimation of the ensemble average (O) of a radiation operator O can 
be obtained by averaging the kernel R[0](x, </>) given in Eq. (|). However, Eq. (|^) needs 
a procedure that exploits the null function equivalence, and is given in Ref. [p8fl . For this 
reason, for simplicity here we use the Richter formula p9| , which expresses the kernels for 
the normally ordered moments as follows 

K.*^*)^^^, (31) 

H n (x) being the Hermite polynomial of order n. We apply the adaptive method to the 
tomographic detection of the most relevant observables: intensity, quadrature and complex 
field amplitude. The optimization method is here particularly useful, as the tomographic 
detection of these observables using the Richter kernel is very noisy p3| , |30| . 
In contrast to the case of matrix elements given in Section [TTlj , here some analytical evalua- 
tions can be carried out. We consider measurements performed on coherent states, squeezed 
vacuum, Fock states and cat superposition of coherent states. It turns out that addition of 
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just few null functions to the Richter kernels generally results in a large improvement of the 
tomographic precision, again with the exception of Fock states where no improvement can 
be obtained. 



A. Intensity 

The tomographic detection of intensity is obtained by averaging the kernel 

R[a!a](x) = 2x 2 - - . (32) 
The vectors b needed for the optimization procedure are given by 

( a t(fc+2)) 



b{ = -R[rta]Fl* = = --^r 1 (33) 

( (at 2 ) _ n 

hi 1 = -R[a^a]F^* = -2x 2 e- i ^+ 1 ) 2 't> - - | ~J~ ™ ~ JJ (34) 

b[ H = -R[tfa]F{ IU = -2x fe W+2 e -*(fcH+2+nH)0 = _ J ^ 1 = . (35) 



1^0 

From Eqs. (§§ and (^) it follows that only F*(x,<j)) and A o 7/ (0) = A o m (0) = exp(2i<p) 
are effective in reducing the variance. We solved analytically the optimization equations 
(PD|) for type-I null functions, and also in this case it turns out that for all the states here 
considered, only the single null function Fq (0) is needed, namely one has 

Ho = b = , V k > 1 . (36) 

The corresponding reduction of variance is easily obtained from Eq. Q2"3|), and is given by 



A 2 [a* a] = ~ (a) 2 ) (a 2 ) . (37) 

Actually, A 2 [a^a] can compensate the leading term of the variance of the original Richter 
kernel which, in turn, is given by 



AR 2 [rfa} = (An 2 ) + - [(a t2 a 2 ) + 2(a f a) + l] . (38) 



This means that the variance of the optimized kernel AA 2 [a^a] becomes much closer to the 
intrinsic intensity fluctuations (An 2 ) than the original noise AR 2 [(tfa\. In order to appreciate 
such noise reduction we compare the two noise ratios 



A A 2 a ta , A A 2 ate, 

Sn R = x J 5n K = J ^-L 1 . 39 

' (An 2 ) V (^ 2 ) 



For coherent states \a) we obtain 
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S "«= \/ 2 + - ( l"l' 2 + — ) »<>!< - \/ 2 +2j^j2 ' ( 40 ) 

that is, from an asymptotically linearly increasing function of \a\ the ratio becomes a con- 
stant Suk — y/2. Similar expressions are obtained for other kind of state: the noise ratio 
saturates to 5rix ~ a/3/2 for either squeezed vacuum and cat states. 

In Fig. ^ results from a Monte Carlo simulation of the tomographic measurement of intensity 
on coherent states show the noise reduction obtained when using the optimized kernel. 
The noise reduction obtained by adding the single null function F (<j)) can be easily evaluated 
also for the generic diagonal moment (a^ n a n ), using the formula 

j 2 *R[<J n a n \(x) = ^-^R[a^ n+1) a n - 1 \(x) , (41) 

which leads to 

6 = -R[rf n a n }e i2 <i> = -^-^(a^V" 1 ) , (42) 

namely A 2 [a^ n a n ] = 2|6 | 2 - We just mention that optimizing the kernel R[a^ 2 a 2 ](x) is useful 
to improve detection of the second order correlation function = {a) 2 a 2 ) / {a) a) 2 . 



B. Quadrature 

The optimization procedure has been tested also on the kernel R{x](x,4>) = 2xcos0, 
corresponding to the measurement of the quadrature operator x = ^(a + al). Similarly to 
the intensity case, the type-II and type-Ill null functions do not play a role in improving 
precision, whereas type-I functions give bk = — 2~ k ~ 1 (a^ 1+k ^) in Eq. (|20|). In this way the 
optimization procedure can be carried analytically also in this case. The results indicate that 
for coherent states it is enough to add the first null function F$ ((/)), whereas for squeezed 
vacuum and cat states only the odd- index functions F 2 / s+1 (a;, 0) contribute to noise reduction. 
In this case the main term is due to Ff(x,<p), whereas higher order functions improve the 
variances only by a few percent. For coherent states the variance reduction from Fq (x, 0) is 
given by 

A 2 [x]= 1 -(a^)(a) = ^\a\ 2 , (43) 
which completely compensates the leading term in the variance of the original Richter kernel 



AR 2 [x] = (Ax 2 ) + i(a f a) + - . (44) 



For squeezed vacuum and cat states the variance reduction due to Fi (x, <f>) is 



A 2 \x] 



2(1- |( G ) |2 + 2<ata» 



|(a)H(at 2 ) + (at 2 ) + i + (ata)) - -|( f r)| 2 



(45) 
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Upon defining the noise ratio Sxk in analogy to Eq.(|39"D 



6x K =J , 46 

V (Ax 2 ) 

from Eqs. (fHf) and (|45|) we get the constant 5xk = V% for coherent states, independently on 
|a| 2 , whereas for squeezed vacuum and cat states the noise ratio saturates to 5xk — a/5/4. In 
Fig. H results from a simulated experiments of tomographic measurement of the quadrature 
on coherent states are shown for | ct | 2 = 3. There the histograms of the original Richter kernel 
and of the optimized kernel are compared. The optimized kernel has a sharper distribution, 
which is peaked at the mean value (x) = For this reason, it is quite obvious that the 
optimized kernel K[x](x,4>) gives a more precise determination of (x) than the the original 
kernel i?[x](a;, <fi). 



C. Field amplitude 



The tomographic kernel for the measurement of the complex field amplitude a is given 
by R[a](x,4>) = 2xe 1 ^, and its fluctuations should be compared with those from the ideal 
measurement of a, which could be achieved by ideal eight-port p2|-p4| or six-port p5|j36 



homodyne detection. The optimization procedure depends on the choice for the definition 
of statistical error for a complex quantity. If one considers the real or the imaginary part 
separately, the procedure coincides with the optimization of the precision in independent 
measurements of two conjugated quadratures. On the other hand, in order to take into 
account both noises jointly, we minimize the quantity 



1 



A,K 2 [a] = - <! \K\n] 



K 2 \ 



(47) 



corresponding to the average of noises for real and imaginary parts, namely the trace of the 
noise covariance matrix. Now, the equivalence class of kernel functions is written as follows 



Af-l Af-l 

K[a](x,(j)) = R[a](x,<p) + (ipF p (x,<p) + z/ p F p *(x,0) 

p=0 p=0 



(48) 



H p and Up being two independent sets of complex coefficients. The optimization procedure 
is similar to the real case, and is reduced to solving the two linear systems 



A/i = b Av = c, 

where c is given by 

c p = -R[6]Fp . 
By inverting Eqs. (f49D , one obtains the noise reduction 

M-l 

A 2 [a] = AT^R - AJO\a\ = £ [b p {A^) gp b* q + c p (A^) pg c* 

p,q=0 



(49) 



(50) 
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Also in the present case it is sufficient to consider only type-I functions. The optimization 
vector b is given by b k = -2- k (a^ 1+ V). Similarly to the case of the quadrature, the 
optimization procedure shows that for coherent states only Fq{4>) is needed, whereas for 
squeezed vacuum and cat states only the odd- index functions F% s+l (x, 0) contribute to noise 
reduction, and the main term comes from F((x,(j)). In this way for coherent states one 
obtains 



A:[a] = -|a| 2 , 

whereas for squeezed vacuum and cat states one has 

1 



Al[a] 



\(a)\ 2 [(a^) + (a^) + - + (a^a))+\(a 2 )\ 2 



(51) 



(52) 



2 (1 - |(a>| 2 + 2(at a >) [ 
Eqs. ( |5TD and (|52|) should be compared with the noise-figure of the original Richter kernel 



AlR{a} = -[2(a^a) + l-\(a)\ 2 } 



and with the intrinsic noise of a generalized measurement of the amplitude 

1 



(A^) = -[(ata) + l-|(a)| 2 ] 



(53) 



(54) 



The noise ratios thus equals Sa^ = 1 for coherent states, whereas saturates to 5a K ~ a/3/2 
for both squeezed vacuum and cat states. Remarkably, for coherent states the heterodyne 
noise is reached, namely tomographic detection has ideal noise. 



V. EFFECTS OF SYSTEMATIC ERRORS 

Throughout this paper the tomographic kernels have been optimized by adding low order 
null functions. Higher order functions oscillate more rapidly. Since the method involves only 
the average of these functions on a small sample of data, fast oscillations in and higher 
power of x would introduce more noise, and including too many null functions would increase 
the error instead of reducing it. In Fig. ||] an example of such pathology is given. 
Another point that should be mentioned is that in the tomographic detection here consid- 
ered the phase is a random parameter in [0, 7r]. A discrete scanning by equally-spaced 
phases would introduce systematic errors [j2T| , |3~lf that would mask the benefits from the 
optimization. Actually, for non-random uniform scanning, the null function -F o (0) has no 
effects when added to phase independent kernels, whereas the other null functions have a 
much reduced effect, and obviously do not eliminate the systematic error due to the finite 
mesh of the deterministic scanning. 



VI. SUMMARY AND CONCLUSIONS 



In this paper we have presented an adaptive method to optimize tomographic kernels, 
improving the precision of the tomographic measurement. The method has been analyzed 
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in detail for coherent states, Fock states, squeezed vacuum, and " Schrodinger-cat" states. 
With the exception of Fock, states the method generally provides a sizeable reduction of 
statistical errors. For coherent states the improvement mainly concerns the small-index 
matrix elements, whereas for squeezed vacuum and cat states also far off-diagonal elements 
are improved. 

The error reduction is much more significant for the measurement of intensity, quadrature 
and field amplitude, where for coherent states, squeezed vacuum, and cat states the ra- 
tio between tomographic noise and uncertainty of the considered observable saturates for 
increasing energy. In this case, we can definitely assert that quantum tomography is a quasi- 
ideal measurement, as it adds only a small amount of noise as compared to ideal detection. 
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FIGURES 
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FIG. 1. Noise reduction in the tomographic measurement of the diagonal elements {n\g\n) of 
the density matrix of a coherent state \a) with intensity \a\ 2 = 5. In (a): noise reduction 7 versus 
the number of added type-I null functions: the full curve represents (0|f?|0), the dashed curve 
and so on, from the top to the bottom. In (b): noise reduction versus the index n of the 
diagonal matrix element for M = 30 added null functions. 
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FIG. 2. Monte Carlo simulation of adaptive tomography of a coherent state with intensity 
|q| 2 = 3. A sample of 5 blocks of 50 homodyne data is used for each of 25 phases (for a total 
number of measurements equal to 6250). The optimization has been performed by adding M = 6 
null functions. In (a) the measured diagonal matrix elements before optimization, and in (b) after 
optimization. The squares indicate theoretical values. 
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FIG. 3. Noise reduction 7 versus indices n and m of the matrix element (m\g\n) for a coherent 
state with intensity \a\ 2 = 5: (a) using only type-I null functions, (b) using only type-II, and (c) 
type-Ill. For all plots M = 10 null functions have been used in the optimization procedure. 
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FIG. 4. Noise reduction for squeezed vacuum with (n) = 4: (a) using only type-I null functions, 
(b) using only type-II, and (c) type-Ill. For all plots M = 10 null functions have been used in the 
optimization procedure. 
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FIG. 5. Noise reduction for the cat-like superposition of coherent states and for the three types 
of null functions in Eq. ( f28|) with a = V3: (a) using type-I null functions, (b) using type-II, and 
(c) type-Ill. For all plots M = 10 null functions have been used in the optimization procedure. 
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FIG. 6. Adaptive tomography of a squeezed vacuum with (n) = 4. The Monte Carlo sample 
includes 5 blocks of 100 data for each of 50 phases (for a total number of measurements equal to 
25000). The optimization has been performed by adding M = 10 type-II null functions, (a) mea- 
sured elements without optimization; (b) with optimization. The squares indicate the theoretical 
values. 
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FIG. 7. Tomographic detection of the intensity on coherent states. The simulated experiment 
has been performed with 15 blocks of 15 data for 15 phases each (for a total number of N = 3375 
measurements). The tomographic result (h) is reported versus the theoretical values \a\ 2 , (a) 
without and (b) with optimization. 
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FIG. 8. Histograms of the kernel functions evaluated on the tomographic outcomes for a co- 
herent state with \a\ 2 = 3. The sample has 50 phases with 100 data each, (a) using the original 
Richter kernel R[x](x,(p); (b) using the optimized kernel K[x](x,(p) . The distribution for the 
optimized kernel is sharper and peaked near the theoretical value (x) = 
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FIG. 9. Monte Carlo simulation of adaptive tomography with a bad choice of the number of 
added null functions. The state under examination is a coherent state with \a\ 2 = 3 and the 
simulated sample of homodyne data contains 5 blocks of 50 data for 25 phases each, for a total 
number of 6250 measurements (as in Fig. Here, the optimization has been performed by adding 
M = 32 null functions. Large fluctuations emerge instead of error reduction. The squares indicate 
theoretical values. 
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TABLES 



TABLE I. Representation table for indices of type- III null functions. 
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